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\ Abstract 

The influence of a frustrated bond on the magnetic properties of a d=3 uniaxial (Ising) b.c.c. 
diluted antiferromagnet, with emphasis in the compound FexZni-xF2, is investigated by a 

i local mean-field numerical simulation. In particular we find that the initial drop of the satura- 

tion staggered magnetization (Ms) with concentration follows a percolation-like phenomenon 

I characterized by an exponent Pp. For the frustrated samples, however, this regime is followed 

by a second one identified by a "long tail" effect such that Ms is zero only at the percolation 

I threshold. Our numerical data also confirms a spin-glass phase near this threshold. 

; PACS. 75.50.L - Spin glasses. 

' PACS. 75.50.E - Antiferromagnetics. 

PACS. 61.43.B - Computer simulations in disordered solids. 
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Diluted antiferromagnets (DAFs) have been largely used as testing ground for the- 
oretical models of random magnets[l]. In spite of this, a complete understanding of 
the microscopic mechanisms behind the rich variety of phases found in these systems 
is still far from being accomplished. The particular nature of the exchange interac- 
tions, the external applied field, as well as the spatial anisotropy give rise to distinct 
magnetic properties in real systems. With respect to the above parameters random 
systems with uniaxial anisotropy and short-range interactions can be classified by 
Ising models belonging to different universality classes. In three dimensions (d=3) 
the insulating DAF compound FexZni-xF2, with a centred tetragonal lattice, is 
one of the most studied experimental realizations of such models[l-4]. Due to a 
strong single- ion anisotropy, pure FeF2 is considered an archetypal Ising system. 
The partial random substitution of Fe"^^ by Zn'^'^ ions in this compound make it to 
behave magnetically as a random- exchange Ising model (REIM). The application of 
a uniform field parallel to its easy direction causes a crossover from REIM to the 
random- field Ising model (RFIM) critical behavior [1,5]. For highly diluted samples, 
random fields of strong magnitudes make the AF long-range order (LRO) unstable 
and nucleates a "glassy" phase in the upper part of the {H,T) phase diagram [6]. In 
the vicinity of the percolation threeshold {xp ^ 0.24) a true Ising spin glass (SG) 
phase was found to occur below a freezing temperature in a sample of 
[4]. For the appearence of a SG behavior randomness and frustration are consid- 
ered essential ingredients . However, the effects of very weak frustrated interactions, 
eventually present in short-ranged random magnets, have been a priori neglected as 
compared with the ones caused by stronger interactions [2]. The opposite limit, i.e., 
systems with weak disorder and strong frustrated interactions are also of consider- 
able interest. 

In a previous paper [3], we presented results of a zero-field (H = 0) numerical 
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simulation which showed the dramatic role played by a weak frustrated interac- 
tion in the antiferromagnetic (Ms) and spin-glass (Q) order parameters of the 
FexZni-xF2, leading the system to the SG phase observed [4] in the strong-dilution 
regime {x « 0.25) of this compound. Here, we cover the influence of a frustrated 
bond in a centred tetragonal lattice of a DAF (e.g. , FexZni-xF2) ai H = 0, 
for an extended range of frustration strenghts in the entire interval of dilution 
{0 < X < 1.0). Ms{x) is shown to be strongly dependent on the magnitude of 
frustration present in the whole range oi x. In particular it follows a percolation- like 
phenomenon characterized by an exponent /3p. For frustrated samples a "long tail" 
effect is identified such that Ms is zero only at Xp. On the other hand, Q{x) is 
almost independent on the frustration strenght, except for x close to Xp where it 
suffers an abrupt increase in magnitude even for the case of a very weak frustration. 
We use a microscopic Hamiltonian suitable to describe a randomly diluted Ising 
antiferromagnet, with short-range exchange interactions Ji, given by: 

= X] Ji^i^i+SeSiSi+Se^ (1) 
<i,6i> 

where Si = ±2 (as in the Fe+^), Ei = 0,1 and i is summed up over the three 
sets of nearest neighbors in the centered tetragonal lattice shown in Fig. 1. We 
take N = 2 • 30^ sites and use periodic boundary conditions. The local mean-field 
(LMF) method consists in solving iteratively the self-consistent equations involving 
the thermally averaged spins rrii obtained through a variational minimization of the 
MF-free energy [7], which gives for Si = 1 

rrii =< >T= 2tanh[(l//c5r) ^ J^^.+j^m.+jJ. (2) 

The usual virtual crystal mean field (VCMF) result is obtained assuming rrii = m 
for every magnetic site. The LMF simulation starts by choosing a random initial 
configuration {Si, Si} in the high-T (paramagnetic) phase. The system is cooled 
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in steps AT = IK through (2) and heated back by the same amount up to a 
temperature of interest, following a convergence criterion: 



tin 
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where n represents the n-th iteration. We measured the staggered magnetization 

MsiT,x) = i2/N)j:£^mi, (4) 

i 

with i summed over a given sublattice, and the spin-glass Edwards- Anderson (EA) 
order-parameter 

Q{T,x) = il/N)j:e,<S]>T, (5) 
j 

with j taken through the whole lattice. Ms and Q are averaged over 50 indepen- 
dent samples. In the simulations of the FexZni-xF2 compound we considered the 
exchange constants Ji, shown in Fig. 1, keeping unaltered the experimental ratios 
ji = J1/J2 = —0.013 and js = J3/J2 = +0.053 measured[8] for FexZni-xF2, such as 
to fix the Neel temperature of the pure system, Tn{x = 1) = 77.8 K. This procedure 
appears to be quite realistic for application in FexZni-xF2i as early pulse experi- 
ments [9] in this system support no a:-dependence of the total exchange interaction 
in the entire concentration range measured (0.2 < x < 1.0). 

The x-dependence of the normalized sublattice magnetization Ms{x)/Ms{x = 1) is 
shown in Fig. 2, varying the frustrated interaction in the whole interval < j3 < 1. 
Data were obtained at T/Tn{x = 1) = 0.18 and H = 0. The full line represents the 
expected result of the VCMF method for the specific case of FeF2, where = 0.053 
. As well known, the VCMF method does not predict the existence of a percolation 
threshold. Moreover, in this method the presence of a weak frustrated interaction 
(js = 0.053) plays no role in the linear dependence of Ms with x. By way of con- 
trast, in the LMF method we observe that even a very weak frustration (ja < 0.053) 
reduces dramatically the values of Ms in the strong dilution regime(x < 0.5), as is 



the case of FexZni-xF2 [3]. However, as can be seen from the J3 = 0.053 curve, a 
weak frustration plays no effect on the AF order parameter Ms of this compound for 
weak dilution (0.6 < x < 1). On the other hand, a strong frustration (jf's 1) may 
change the nature of the long range order even in the pure system, x = 1. Ms{x) 
vanishes at the percolation concentration Xp = 0.247 ± 0.010, determined by the 
projection to the x-axis of the best fit curve from the = data. This is in very 
good agreement with the value Xp = 0.243 ± 0.010, obtained[10] by series methods 
for non- frustrated b.c.c. lattices. In the LMF method Ms becomes smaller than the 
VCMF values for x below a frustration-dependent value x{js)- We perceive from 
the data of Fig. 2 that x{j^) increases for increasing ja, but reaches the saturation 
value a; = 1.0 for j3 < 1.0 . For non-zero frustration strenghts the Ms curves display 
tails toward the percolation threshold, after sudden initial drops starting at ^(ja). 
Close to the percolation threshold, the antiferromagnetic order parameter have the 
asymptotic form Ms ~ (x — Xp^p, where /3p is the percolation exponent associated 
to the order parameter. The best fitting of this law, with = 0, gives the value 
/3p = 0.63 ±0.11 . We notice the same behavior in the frustrated systems for which 
we find (3p = 0.59 ± 0.10 for js = 0.053 (FeFz), and Pp = 0.62 ± 0.09 for h = 0.265. 
Series expansion[ll] and renormalization group[12] methods yield values for /3p be- 
tween 0.4 to 0.5 . One should stress, however, that Ms vanishes only ai x = Xp for 
the whole interval < js < 1.0 . 

The normalized Edwards- Anderson order parameter Q{x)/Q{x = 1) is plotted ver- 
sus x, for several values of ja in Fig. 3. The linear dependence of Q{x) with x, found 
for X > 0.4, is not disturbed even for the strongest magnitude of frustration used in 
this work (ja = 1-0). However, close to the percolation threshold even an evanescent 
frustration, as is the case of FexZni-xF^^ causes a pronounced increasing in Q{x)^ 
as shown in the inset of Fig. 3. The values Ms[x = 0.25) « and Q[x = 0.25) 7^ 
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found in this LMF simulations for j-^ = 0.053 explains [3] the spin glass phase exper- 
imentally detected in For x < 0.24 however, both Ms and Q values 
are negligibly small (see Figs. 2 and 3, respectively) for the whole T interval mea- 
sured in the present work (2K < T < 150K). So, the LMF method does not support 
a true spin-glass phase {Q 0, Ms ~ 0) for very low concentrations of magnetic 
ions, in conformity with the short-range character of the exchange interactions. 
In the LMF approach Tn{x) is determined by the inflection point of the Ms{x) 
versus T curve (see Ref. [3]). The x dependence of the normalized Neel tempera- 
ture Tn[x)/Tn{x = 1) , obtained by VCMF and LMF methods, is compared with 
experimental results for the FexZni-xF2 compound in the phase diagram of Fig. 4. 
The LMF predictions for Tn{x) are in good agreement with experimental data in 
the compound FexZni-xF2 for the whole range of x, in spite of the increasing er- 
ror bars in the vicinity of the singular point Xp. This large uncertainty is due to 
numerical difficulties to obtain Ms at very low temperatures {T/Tn{x = 1) < 0.17) 
in the strong-dilution regime (0.24 < x < 0.31). Notice that Tn(x) is virtually 
unchanged even by the presence of moderate strenghts of frustration (see Fig. 2). 
This is explained by the fact that in spite of a sudden reduction on the magni- 
tude of the antiferromagnetic order parameter Ms{x) when a frustrated interaction 
(0 < is < 1-0) is present, Ms{x) ^ for all x> x^ (see Fig. 2). However, we cannot 
exclude the possibility that the inclusion of thermal fluctuations may destroy such 
tiny residual Msix)^ leading Tn{x) to vanish for x > Xp, as experimentally observed 
in FexZni-xF2 [13]. 

Our reported results show that the presence of a small frustrated interaction plays 
no role on the magnetic properties of a DAF system for weak dilution. However in 
the moderately and strongly diluted regimes it causes major effects on a variety of 
magnetic properties, some of which have been identifled in this work. 
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Figure Captions 

Fig.l b.c.c. magnetic structure, illustrating the first three nearest-neighbors exchange 
interactions between magnetic ions. |^ 

Fig. 2 x-dependence of normalized sublattice magnetization M5(x), atr/r/v(x = 1) = 
0.18 and H = 0, for various exchange ratios = J3/ J2, measured by LMF 
simulations (symbols) and VCMF method (full line). 

Fig. 3 x-dependence of normalized spin-glass Edwards- Anderson order parameter Q(x), 
at T/T]\i{x = 1) = 0.18 and i7 = 0, for various exchange ratios ja = 
measured by LMF simulations. 

Fig. 4 Tn{x)/Tn{1) vs. x using VCMF (full line) and LMF simulations in the frus- 
trated (•) and non-frustrated (□) cases, H = 0. Experimental data (A) of 
FexZni-xF2 are plotted for comparison (Ref. [13]). Xp = 0.247 ± 0.010 indi- 
cates the percolation threshold. 



^Unfortunately we don't have a ps file of this figure. However, it is also reproduced in Ref. [3] 
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